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Abstract 

o : 

. The formation of K~ and K° condensation in /3-equilibrated hyperonic mat- 

£H ! ter is investigated within a relativistic mean field model. In this model, 

baryon-baryon and (anti)kaon-baryon interactions are mediated by the ex- 
change of mesons. It is found that antikaon condensation is not only sensitive 
■ to the equation of state but also to antikaon optical potential depth. For large 

values of antikaon optical potential depth, K~ condensation sets in before the 
appearance of negatively charged hyperons. We treat K~ condensation as a 
\ first order phase transition. The Gibbs criteria and global charge conserva- 

tion laws are used to describe the mixed phase. Nucleons and A hyperons 
\& ■ behave dynamically in the mixed phase. A second order phase transition to 

K° condensation occurs in the pure K~ condensed phase. Along with K~ 
, condensation, K° condensation makes the equation of state softer thus re- 

^ (— i suiting in smaller maximum mass stars compared with the case without any 

Q-f condensate. This equation of state also leads to a stable sequence of com- 

pact stars called the third family branch, beyond the neutron star branch. 
> ■ The compact stars in the third family branch have different compositions and 

smaller radii than that of the neutron star branch. 
PACS: 26.60.+C, 21.65.+f, 97.60.Jd, 95.30.Cq 

•i-H . 

I. INTRODUCTION 

There is a growing interest to understand whether a stable sequence of compact stars 
could exist in nature beyond a neutron star branch. ^From observations, two families of 
compact stars - white dwarfs and neutron stars, are known to us. The physical reasons for 
stability in white dwarfs and neutron stars are different. In white dwarfs, it is the Fermi 
pressure of degenerate electrons that stabilizes the stars. After the white dwarf branch, 
there exists no stable star having central densities in the range 10 9 — 10 l4 gm/cm 3 as is 
evident from the mass-radius relationship obtained by the Tolman-Oppenheimer-Volkoff 
(TOV) equations Q. The stability is regained in the neutron star branch. In this case, 
the Fermi pressure of degenerate and interacting baryons is responsible for the stability of 
neutron stars. The question is what is next to a neutron star branch - again an unstable 
region followed by a stable configurations of super dense stars! If this picture is true, what 
could be the mechanism behind the stability of a possible third family of compact stars? 
Gerlach [0] first pointed out that a third family of compact stars could be a possibility. He 
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noted that a jump in the equation of state (EoS) or equivalently a large discontinuity in 
the speed of sound beyond the central density corresponding to the maximum neutron star 
mass might result in a stable family of super dense stars beyond the neutron star branch. 
However, Wheeler et al. showed that there was no stable compact stars beyond a neutron 
star branch for a smooth EoS [[|. 

Neutron stars are very useful laboratories to investigate the properties of dense matter 
At the core of a neutron star, the matter density could exceed by a few times normal 
nuclear matter density. Various exotic forms of matter such as hyperonic matter, quark- 
hadron mixed phase and Bose-Einstein condensation of strange particles may appear there. 
Earlier it was found that there might be kinks in the EoS because of the appearance of an 
exotic form of matter in the high density regime. Consequently there is a discontinuity in 
the speed of sound. This motivated various groups to explore whether physical processes like 
a first order phase transition from hadronic matter to quark matter or hadronic matter to 
hyperonic matter could give rise to an EoS leading to a third family of compact stars It 
was shown by various authors that a first order quark-hadron phase transition could produce 
an EoS such that the neutron star branch is terminated by the softening in the EoS due to the 
mixed phase and the third family of compact stars attained a pure quark phase in the core 
. Similarly, a first order phase transition from hadronic matter to hyperonic matter led 
to stable solutions beyond the neutron star branch for certain parameters sets of hyperon- 
hyperon interactions 0. It was demonstrated in those calculations that non-identical stars 
of same mass but distinctly different radii and compositions could exist because of partial 
overlapping mass regions of the neutron star branch and the third family branch. These 
pairs are known as "neutron star twins" 

Besides quark-hadron phase transition and the formation of hyperonic matter, Bose- 
Einstein condensation of antikaons may occur in the /3-stable and charge neutral dense 
matter. This antikaon condensation could be a first or second order phase transition || . It 
was first demonstrated by Kaplan and Nelson [|TUJ using a chiral SU(3) l x SU(3)r Lagrangian 
that i^~meson may undergo Bose-Einstein condensation in dense hadronic matter formed in 
heavy ion collisions. The strongly attractive .ff "-nucleon interaction lowers the effective mass 
of K~ mesons in dense matter. Consequently, the in-medium energy of K~ mesons decreases 
and s-wave condensation sets in when the in-medium energy of K~ mesons equals to its 
chemical potential. Later, this chiral model was applied to the study of K~ condensation in 
the core of neutron stars |TT|-Q. Also, the traditional meson exchange picture where baryons 
and (anti)kaons interact through meson exchange, was extensively exploited to investigate 
antikaon condensation in dense matter relevant to (proto)neutron stars HJl5|-pO|. The onset 
of K~ condensation in neutron star matter depends on the nuclear equation of state and 
also on the depth of antikaon optical potential. 

The threshold density of K~ condensation in nucleons-only star matter is 3-4 times 
normal nuclear matter density. The net effect of K~ mesons in the condensate in neutron star 
matter is to maintain charge neutrality replacing electrons and soften the EoS resulting in the 
reduction of maximum mass of the star |^|,|T2|,|rS|-0,|T^,^] . In the presence of hyperons, K~ 
condensation was delayed to higher density. Recently, we have investigated the formation of 
K° condensation along with K~ condensation in dense nuclear and hyperonic matter relevant 
to (proto)neutron stars within Relativistic Mean Field (RMF) models fL9| , 20| . It was found 



that the onset of K condensation always occurred later than that of K condensation. 
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With the onset of K° condensation, abundances of neutrons and protons become equal in 
the high density matter [|l^,^(| • I n the presence of hyperons, the thresholds of K~ and K° 
condensation are shifted to higher densities. Along with K~ condensate, K° condensate 
makes the EoS much softer resulting in smaller maximum mass star compared to the case 
without any condensate. 

In this paper we investigate the effects of antikaon condensation on the equation of state 
within a RMF model where baryons and (anti)kaons are interacting via meson exchanges. 
Also, we study the properties of the third family of compact stars with this EoS. The paper 
is organised in the following way. In section 2, we describe the RMF model of strong 
interactions and different phases of dense matter. In section 3, the parameters of the model 
are discussed. Results of our calculation and their relations to the third family of compact 
stars are explained in section 4. Section 5 deals with summary and conclusions. 



II. FORMALISM 

In this calculation, we describe a first order phase transition from hadronic matter to the 
antikaon condensed phase in compact stars. We adopt a relativistic field theoretical model 
to describe the pure hadronic matter, pure antikaon condensed matter and the mixed phase. 
The constituents of matter - n,p, A, S + , £~, E°, 5~, H° of the baryon octet and electrons 
and muons, in the hadronic phase are to satisfy beta-equilibrium and local charge neutrality 
condition. The strong interaction between baryons is mediated by the exchange of scalar 
and vector mesons. The model is also extended to include hyperon- hyperon interaction 
through two additional hidden-strangeness mesons- scalar meson / (975) (denoted hereafter 
as a*) and the vector meson 0(1020) ||17|| . Therefore the Lagrangian density for the pure 
hadronic phase is given by 

£b = J2^b (hfid^ -m B + g aB cr - g^Bl^ - g pB l^B • P M ) ipB 

B 

+ l - (d^Pa - mla 2 ) - U(a) 

--ov^ + \mlu^ - ^ ■ fT + l -m 2 p p^ ■ p» + Cyy . (1) 

Here ips denotes the Dirac spinor for baryon B with vacuum mass m# and isospin 
operator t#- The scalar self-interaction term |2T| 



is 



U{&) = \92^ + . (2) 

The Lagrangian density (Cyy) responsible for hyperon- hyperon interaction is given by, 

Cyy = ^24>b (ga*B<?* - g^Bi^) ^b 

B 

+~ (d»a*d»a* - ml,a* 2 ) - \<j>^ + \m\^ ■ (3) 

At the interior of compact stars, the generalised /3-decay processes may be written in 
the form B\ — > B 2 + / + V\ and B 2 + / — > B\ + V[ where B\ and B 2 are baryons and 1 
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is a lepton. These weak processes conserve baryon number and electric charge and are in 
chemical equilibrium. Therefore the generic equation for chemical equilibrium condition is 

tk = bifin ~ We , (4) 

where fi n , \x e and /ij are respectively the chemical potentials of neutrons, electrons and i- 
th baryon with = (k F . + m* 2 ) 1/2 + g^ujQ + g^cpo + hi9piPoz and 6, and are baryon 
and electric charge of ith baryon respectively. The above equation implies that there are 
two independent chemical potentials p n and /x e corresponding to two conserved charges i.e. 
baryon number and electric charge. In neutron stars, electrons are converted to muons by 
e~ — > /i~ + -\- v e when the electron chemical potential becomes equal to the muon mass. 
Therefore, we have /i e = /i M in a neutron star. In the pure hadronic phase, the total charge 
density is 

Q h = 9Bn% -n e -n^ = (5) 

B 

where is the number density of baryon B in the pure hadronic phase and n e and n M are 
charge densities of electrons and muons respectively. 



Solving the equations of motion in the mean field approximation |22| along with effective 
baryon masses (m*) and equilibrium conditions we immediately compute the equation of 
state in the pure hadronic phase. The energy density (e ) is related to the pressure {P h ) in 
this phase through the Gibbs-Duhem relation 

P h = J2^-e h - (6) 

i 

Here \L{ and rij are chemical potential and number density for i-th species. 

The pure antikaon condensed phase is composed of baryons, leptons and antikaons. In 
this phase, the constituents of matter are in chemical equilibrium under weak interactions 
and maintain local charge neutrality. The baryon-baryon interaction in the pure condensed 
phase is described by the Lagrangian density as given by Eq. (1). Earlier it was shown that 
nucleons in the pure nuclear and antikaon condensed matter behaved differently because 
of their dynamical nature [[J. It was attributed to different mean fields which nucleons 
experienced in the pure phases. We adopt a relativistic field theoretical approach for the 
description of (anti)kaon-baryon interaction ||. Here, baryon-baryon and (anti)kaon-baryon 
interactions are treated on the same footing. We also include K° condensation in this 
calculation and it occurs as a second order phase transition in the pure K~ condensed 
phase. The Lagrangian density for (anti)kaons in the minimal coupling scheme is, 

C K = D^KD^K - m*£KK , (7) 

where the covariant derivative = d /Ji +iguK UJ ij, + ig ( f ) K<f>iJ. J rWpK^K 4 Pu- The isospin doublet 
for kaons is denoted by K = (K + , K°) and that for antikaons is K = (K~ , K°). The effective 
mass of (anti)kaons in this minimal coupling scheme is given by 

m * K = m K - g aK a - g a * K o* , (8) 



where is the bare kaon mass. In the mean field approximation [22|| adopted here, the 



meson fields are replaced by their expectation values. The mean meson fields are denoted 
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by er, a*, ooo, 0o and po3- The dispersion relation representing the in-medium energies of 
K = (K~ , K°) for s-wave (k = 0) condensation is given by 

v k -,k° = m K~ 9u>ku - fiv>x0o =F -z9 P Kpm , (9) 

where the isospin projection J 3 ^ = =Fl/2 for the mesons K~ (— sign) and K° (+ sign) 
are explicitly written in the expression. Since the a and u fields generally increase with 
density and both the terms containing a and u fields in Eq.(9) are attractive for antikaons, 
the in-medium energies of K mesons decrease in nuclear medium. On the other hand, in 
nucleon-only matter p 3 = n p — n n (n p and n n are the proton and neutron densities) is 
negative; thus the p-meson field favors the formation of K° condensation over that of K~ 
condensation. In hyperonic matter, the repulsive <ft nieson term may delay the onsets of 
antikaon condensation ||17|| . The in- medium energies of kaons K = (K + , K°) are given by, 

uk+, ko = m* K + g^ K u + g^K^o ± ^g P KPo3 ■ (10) 

It is evident here that kaon condensation may be impossible in /3-equilibrated matter because 
the a;-meson term is repulsive for kaons and dominates over the attractive cr-meson term at 
higher densities. However, the attractive meson term may decrease kaon energies in the 
presence of hyperons |I7| . 

The meson field equations in the pure condensed phase are derived from Eqs. (l)-(3) 
and (7) as 

m 2 a a = + J29*Bn B : ' s + g (7 Kj2 n K > ( n ) 

B K 

m f a * = J2g**Bn B : ' S + g a * K J2 n K , ( 12 ) 
b K 

m ^o = Y g^sriB - g^Kj2 n K , (13) 

B R 

m l<l>0 = Y 9^BUb -g$Kj2 n K > ( 14 ) 

B K 

m pP03 = J2 QpBhBTlB + g pK J2 hRriK ■ (15) 
b k 

Here the scalar and number density of baryon B in the antikaon condensed phase are re- 
spectively 

R c 2J B + 1 r*F B m* B 2 



Hb ' ~^2^k W+^W k dk} (16) 

4 = (2J B + l)^f, (17) 

with effective baryonic mass m* B = tub — g a BO ' — g<r* b&* ', Fermi momentum kp B , spin Jb, and 
isospin projection I 3B - Note that for s-wave K condensation, the scalar and vector densities 
of antikaons are same and those are given by M 
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n K -, Ro = 2 (cojc-jta + g U KWo + 9^K<h ± ^BpKfHaj KK = 2m* K KK . (18) 

The total energy density in the antikaon condensed matter has contributions from baryons, 
leptons and antikaons and is given by, 

K 1 22,1 3,1 4,1 2 *2,1 22,1 2 i2 , ^ 22 

+ E^4^ f kFB (k? + m$) 1 /*k? dk + Y,- 2 [ KFl (k 2 + mWk*dk + e 3 , (19) 

^ /7T^ JO j 7T JO 

where I goes over electrons and muons and the energy density for antikaons is 

e 3 = m* K (n K - + n^o) . (20) 

Since antikaons form s-wave Bose condensates, they do not directly contribute to the pressure 
so that the pressure is due to baryons and leptons only 

r>R ^ 2 2 1 3 1 4 1 2 *2 1 1 _2 2 1 1 _2 j.2 1 1 2 2 

P =-^ m o° ~^92<y --9z<? -^m a ,a + -m w u + -m^ + -m p p 03 
1 ^ 2J B + 1 rkp B fc 4 dk i v — [ KFl fc4 dk 

The effect of antikaons in the pressure term is through the meson fields, which change due 
to the presence of additional antikaon source terms in the equations of motions in Eqs. (11)- 
(15). In the absence of those source terms for antikaons, we retain the equations of motion 
for meson fields in the pure hadronic phase. 

With the onset of K condensation, other strangeness changing processes may occur such 
as, N # N + K and e" # K~ +u e , where N = (n,p) and K = (K~,K°) denote the isospin 
doublets for nucleons and antikaons, respectively. The requirement of chemical equilibrium 
yields 

V-n ~ A*p = 1*>k- — f^e 1 (22) 
(*k° = 0, (23) 

where fix- and //g-o are respectively the chemical potentials of K~ and K°. The above 
conditions determine the onsets of antikaon condensations. When the effective energy of 
K~ meson (cuj^-) equals to its chemical potential (/z^-) which, in turn, is equal to fi e , a 
K~ condensate is formed. Similarly, K° condensation is formed when its in-medium energy 
satisfies the condition u^o = /i^o = 0. 

The total charge density in the antikaon condensed phase is 

Q K = Yl VbUb ~ n K - n e - n M = (24) 

B 

Now we describe the mixed phase of hadronic matter and K~ condensed matter. In 
earlier works, the Maxwell construction was employed for a first order phase transition 



to antikaon condensation [|i2| , p!6| , p!9| , p0| . It was argued that the Maxwell construction was 
adequate for a system with only one conserved charge PJ23[|. However, neutron star matter 
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has more than one conserved charge. As evident from the generic equation of /3-equilibrium 
(Eq. (4)) there are two independent chemical potentials namely fi n and \x e connected to 
the conservation of baryon number and electric charge, respectively. Therefore, a Maxwell 
construction in this case will result in a discontinuity in one of the chemical potentials 
PJ23|j. It was shown by Glendenning |],^3|] that the Gibbs conditions for thermodynamic 
equilibrium along with global conservation laws would be required to describe the mixed 
phase. For compact star matter the Gibbs phase rules read, 

P h = P R , (25) 

A = /*§■ (26) 

where jj^ and /z^ are chemical potentials of baryon B in the pure hadronic and K~ con- 
densed phase, respectively. The conditions of global charge neutrality and baryon number 
conservation are imposed through the relations 

(i - x)Q h + xQ R = 0, (27) 

n B = (1 - x)n h B + xn% , (28) 

where % is the volume fraction of K~ condensed phase in the mixed phase. The total energy 
density in the mixed phase is 

e=(l- x )e h + X e R . (29) 



III. PARAMETERS 

In the effective field theoretic approach discussed here, knowledge of three distinct sets of 
coupling constants for nucleons, kaons and hyperons associated with the exchange of scalar 
a, isoscalar-vector lu, and vector-isovector p mesons and two additional hidden-strangeness 
mesons, scalar meson (a*) and vector meson <p are required. The nucleon- meson coupling 
constants generated by reproducing the nuclear matter saturation properties are taken from 
Glendenning and Moszkowski of Ref. p4|. This set is referred to as GM1 and its parameters 



are also listed in Table I. 

The vector coupling constants for hyperons are determined from SU(6) symmetry as, 



1 1 1 

1 

= 9p= = g P N', 9 P A = 0, 



2V2 , , 

2fiV,A = 25y, s = = —g^N- (30) 

The scalar meson (cr) coupling to hyperons is obtained from the potential depth of a hyperon 
(Y) in the saturated nuclear matter 

U Y ( U o) = -g<rY<? + guYUQ. (31) 
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The analysis of energy levels in A-hypernuclei suggests a well depth of A in symmetric 
nuclear matter U^(no) = —30 MeV [25,26]- On the other hand, recent analysis of a few 



H-hypernuclei events in emulsion experiments predicts a H well depth of US (no) = —18 MeV 
p7| , |28[| in normal nuclear matter. However, the situation for the £ potential depth in normal 
nuclear matter is very unclear. The only known bound S-hypernucleus is %He The 



most updated analysis of £~ atomic data indicates a strong isoscalar repulsion in E-nuclear 
matter interaction J3UJ. Therefore, we use a repulsive £ well depth of U^(uq) = 30 MeV 
30| in this calculation. 

The cr*-Y coupling constants are obtained by fitting them to a well depth, Uy (no), for 
a hyperon (Y) in a hyperon (Y ) matter at nuclear saturation density |l^,^TJ. It is given as 

UP (n ) = (n ) = 2C/^(n ) = 2U ( A A \n ) = -40 (32) 

Nucleons do not couple to the strange mesons i.e. g a *N = g^N = 0. 

Now we determine kaon-meson coupling constants. According to the quark model and 
isospin counting rule, the vector coupling constants are given by 

QuK = ~q9uN and g pK = g P N ■ (33) 

The scalar coupling constant is obtained from the real part of K~ optical potential depth 
at normal nuclear matter density 

Ur (no) = ~9aK<y - g^K^o ■ (34) 



The strange meson fields also couple with (anti)kaons. Following Ref.[ ||17|| 1, the cr*-K 
coupling constant is determined from the decay of / (925) as g a *K = 2.65, whereas the vector 
(f> meson coupling with (anti)kaons is obtained from SU(3) relation as \[2g^ K = 6.04. 

It was observed in previous calculations that antikaons experienced an attractive poten- 



tial and kaons had a repulsive interaction in nuclear matter |30|j34 -p6l . The analysis of K 



atomic data in the hybrid model [32| yielded the real part of the antikaon optical potential 



as large as Ug = — 180±20MeV at normal nuclear matter density but it was repulsive at low 
density in accordance with the low density theorem. Also, the antikaon potential depth in 
the coupled channel calculation |33| was found to be Ur = —100 MeV whereas the chirally 
motivated coupled channel approach [33J] gave rise to a potential depth of Ur = —120 MeV. 
The different treatments of A (1405) resonance which is considered to be an unstable KN 
bound state just below K~p threshold, may be responsible for the wide range of antikaon 
optical potential depth values in various calculations. In Table II, we list kaon-cx meson 
coupling constant g a K, for a set of antikaon optical potential depths starting from -100 MeV 
to -180 MeV. 



IV. RESULTS AND DISCUSSION 



The conversion of nucleons to hyperons is energetically favorable in (3- equilibrated and 
charge neutral dense matter. Hyperons first appear in dense matter around (2-3)no- Neg- 
atively charged hyperons quench the growth of electron chemical potential. Also the equa- 
tion of state is softened in the presence of hyperons. Therefore, it was shown in various 
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calculations neglecting hyperon-hyperon interaction that the formation of K condensate 



was postponed to higher densities il3,16|, 19| , p0 |. As the matter is hyperon-rich in the high 



density regime, hyperon-hyperon interaction becomes important. This interaction may be 
accounted by including two additional mesons a* and <fi. These strange mesons also couple 
with (anti)kaons. It follows from strangeness numbers of hyperons and K mesons that the 
<p field is repulsive for hyperons as well as for antikaons. The additional attraction due to a* 
field makes the EoS softer. On the other hand, the repulsive contribution of <ft field becomes 
dominant at higher densities and a stiffer EoS is obtained. For calculations with GM1 set 
and various values of |t^-( n o)| < 160 MeV, we find that the threshold of K~ condensation 
in /3-equilibrated hyperonic matter is shifted to very high density (7.5n ) where effective nu- 
cleon mass becomes small [|19|,[2(| • This is attributed to the early appearance of S~ hyperons 
which diminish the electron chemical potential and the presence of repulsive <fi field in the 
expression of energy for K~ and K° mesons (Eq.(9)). 

Now we perform our calculation using GM1 set and antikaon optical potential depth of 
Uj^(n ) = —160 MeV. In Figure 1 we exhibit the particle fraction of /3-equilibrated matter 
containing baryons, electrons, muons and K~ mesons. In the pure hadronic phase where 
local charge neutrality condition is imposed, abundances of nucleons, electrons and muons 
increase with density. Here, charge neutrality is maintained among protons, electrons and 
muons. With the onset of K~ condensation, the mixed phase begins at 2.23n . This mixed 
phase is determined by the Gibbs phase rules for pure hadronic and antikaon condensed phase 
in thermodynamic equilibrium and by global baryon number and electric charge conservation 
laws. We find that A hyperon is the first strange baryon to appear in the mixed phase at 
2.51no- The total baryon density in the mixed phase is the sum of two contributions from 
the hadronic and antikaon condensed phase weighted with appropriate volume fractions. As 
soon as K~ condensate is formed, it rapidly grows with density and replaces electrons and 
muons. Being bosons, K~ mesons in the lowest energy state are energetically more favorable 
to maintain charge neutrality than any other negatively charged particles. Consequently, 
the proton density becomes equal to the density of K~ condensate. Also, the density of 
A hyperon increases with density in the mixed phase. On the other hand, the neutron 
density decreases in the mixed phase. The reason behind it may be the creation of more 
protons in the presence of K~ condensate and also the growth of A hyperons at the expense 
of neutrons. The mixed phase terminates at 4.0no. Heavier hyperons appear in the pure 
condensed phase. As S~ hyperons start populating this phase at 6.77no, the density of K~ 
condensate falls. However, the rising behavior of K~ condensate is regained with the onset 
of S + hyperons at 7.86n . It is to be noted that S~ hyperons do not appear in this case. 

Besides K~ condensation, we now investigate the role of K° condensation in (3- 
equilibrated and charge neutral matter consisting of nucleons, hyperons, electrons and 
muons. The abundances of various species are shown in Figure 2. The composition of 
the pure hadronic and mixed phase in this case is similar to that of Fig. 1. Once the mixed 
phase is over, K° condensation occurs at ~4.1n . We consider K° condensation as a second 
order phase transition. The antikaon condensed phase including K° condensate along with 
K~ condensate is different from that of Fig. 1. With the appearance of K° condensate, 
neutron and proton abundances become equal and a symmetric matter of nucleons is formed 
in this phase. The density of K° condensate increases with baryon density uninterruptedly 
and even becomes larger than the density of K~ condensate. At higher densities, more 
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hyperons start populating the system. As soon as negatively charged hyperons - S~ and 
£~ appear, the density of K~ condensate is observed to fall drastically. This is quite ex- 
pected because it is energetically favorable for particles carrying conserved baryon numbers 
to achieve charge neutrality in the system. Leptons or mesons are no longer required for 
this sole purpose. Moreover, lepton number or meson number is not conserved in the star. 
The system is dominated by K° condensate in the high density regime. 

One important aspect of this calculation is the dynamic behavior of baryons in pure 
hadronic and antikaon condensed phase. Earlier, this dynamical behavior of nucleons was 
noted by Glendenning and Schaffner ||. Besides nucleons, we find similar properties of A 
hyperons which appear in the mixed phase. It is evident from Eq. (11) and Eq. (12) that 
baryons have different effective masses in the hadronic and condensed phase due to different 
a and a* fields. In the pure hadronic phase, the effective mass of a baryon decreases with 
density. With the onset of the mixed phase, a second solution for the effective baryon mass 
emerges. This is the effective baryon mass in the condensed phase. The value of the effective 
baryon mass in the pure hadronic phase is always higher than that of the antikaon condensed 
phase. In the condensed phase, the effective mass of a baryon also decreases with density. 

Equation of state, pressure(P) versus energy density (e), for /3-equilibrated and charge 
neutral matter with different compositions are displayed in Figure 3. Here the long-dashed 
curve stands for hyperonic matter without any antikaon condensate. The dotted curve 
represents hyperonic matter including K~ condensate whereas the solid curve corresponds 
to hyperonic matter with both K~ and K° condensate. The EoS denoted by the dashed-dot 
curve is same as given by the solid except no H hyperon is present in the system in this 
case. Two kinks in the solid curve marks the beginning and end of the mixed phase where 
pure hadronic and K~ condensed phase are in thermodynamic equilibrium as dictated by 
the Gibbs conditions and global conservation laws. Immediately after the termination of 
the mixed phase at 4.0no, K° condensation sets in at ~ A.Iuq. We have treated here K° 
condensation as a second order phase transition. Also we find that the pressure increases 
continuously with energy (solid curve) after the formation of K° condensation. Earlier it was 
noted that the antikaon condensation could be a second order phase transition depending 
on the strength of antikaon optical potential depth and parameter sets of different models 
H Jl9|j20[| . Though the dotted curve has overlap with the solid curve up to the mixed phase, 
it becomes stiffer in the high density regime. On the other hand, the dashed-dot curve 
becomes stiffer compared to the solid curve because there is lesser number of degrees of 
freedom with the absence of H hyperons in the former case. With the appearance of both 
K~ and K° condensation in hyperonic matter, the overall EoS (solid curve) is softened due 
to the strong attraction imparted by antikaon condensation. 

We have used the results of Baym, Pethick and Sutherland |37| to describe the crust of 
a compact star composed of leptons and nuclei for the low density (jib < 0.001/m~ 3 ) EoS. 
In the mid density regime (0.001 < ub < 0.08/m~ 3 ) the results of Negele and Vautherin 
38| are taken into account. Above this density, an EoS calculated in the relativistic model 



has been adopted. 

Now we present the results of static structures of spherically symmetric compact stars 
calculated using TOV equations. The static compact star sequences representing the stellar 
masses M/M & and the corresponding central energy densities e c are shown in Figure 4. 
We denote the sequence of neutron stars with K~ condensation in hyperonic matter by 
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the dotted curve. Here, the maximum neutron star mass is 1.649M occurring at central 
density 7.7n Q . The maximum neutron star mass consisting only of hyperons is 1.789M 
corresponding to the central density 5.16no (not shown in the Figure). The solid curve 
represents compact stars with both K~ and K° condensation in hyperonic matter. When 
K° condensation along with K~ condensation is taken into account, the limiting mass of 
the neutron star branch is attained at much earlier central density 4.49no- In fact, the 
limiting mass star appears very close to the upper boundary of the mixed phase and has a 
smaller value 1.571M Q than the situation with only K~ condensation. This reduction in the 
maximum mass of the star may be attributed to the softening of the EoS in the presence 
of K~ and K condensate. After the stable neutron star branch, there is an unstable region 
followed by a sequence of compact stars in Fig. 4. This branch of super dense stars beyond 
the neutron star branch is called a third family of compact stars There are two third 

family branches in Fig. 4. The lower curve depicts the third family of compact stars with 
both antikaon condensates in hyperonic matter (case I) whereas the upper curve represents 
the third family of compact stars including both condensates in hyperonic matter except 
H hyperons (case II). The maximum masses of compact stars in the third family branches 
denoted by the lower curve and upper curve are 1.552M and 1.561M corresponding to the 
central densities 8.09n and 9.20n respectively. It is the behavior of the EoS including both 
K~ and K° condensate in the mid and high density regime that results in a third family of 
super dense star. The possibility of a third family of compact stars was already discussed in 
the context of quark-hadron phase transition |||| and hadronic matter to hyperonic matter 
phase transition 0. 

We have found a new sequence of compact stars called the third family of compact stars 
beyond the neutron star branch as is evident in Fig. 4. So far, we have discussed the 
hydrostatic stability of the stellar sequence in the third family using TOV equations [|TJ. 
In Fig. 4, we note that the neutron star branch or the third family branch has positive 
slope i.e. dM/de c > 0. Though the positive slope of a steller sequence is a necessary 
condition for stability, dynamical stability requires the analysis of the fundamental mode 
of radial vibration []5|,|6,39- f4l|l . In the third family branches in Fig. 4, we find that each 
configuration is stable because the squared frequency of the fundamental mode associated 
with it is positive. Similar results were obtained in the study of compact stars beyond the 
neutron star branch using quark-hadron phase transition by various groups [f^,f|. 

We exhibit the mass-radius relationship for different EoS in Figure 5. The maximum 
mass star composed of hyperons and no condensate (not shown here) has a radius of 12.8 
Km whereas it is 10.9 Km for the maximum mass neutron star including hyperons and K~ 
condensate (dotted curve). Two cases where we obtain stable solutions beyond the neutron 
star branch are shown in the inset. The solid curve denotes the sequences of compact stars 
having hyperons and K~ and K° condensate. In the neutron star branch, the maximum 
mass star has a radius of 12.8 Km. For the third family branches, the lower curve corresponds 
to case I and the upper curve denotes case II. The radius of the maximum mass super dense 
star on the lower curve is 10.7 Km whereas that of case II is 10.2 Km. In this calculation, 
we find that the existence of non-identical stars of same mass may be possible because of 
partial overlapping of the neutron star and the third family branch. These pairs are known 
as "neutron star twins" |5]||. Also, we observe that the high density twin is a very compact 
object with a different composition and smaller radius than its low density counter part. It 
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is to be noted that the radius of the maximum mass star with K~ condensate in hyperonic 
matter is comparable to that of the third family considered here. This implies that only 
one star with small radius cannot ensure the existence of a third family of compact stars. 
Schertler et al. || argued that two stars - one in the neutron star branch and the other in 
the third family branch with similar masses but with a radius difference of a few kilometer 
might reveal the existence of a third family of compact stars. In this case, the mass difference 
(AM) of two stars is to be much less than the radius difference (AR) i.e. AM << AR. In 
our calculation, we note that AR is almost two orders of magnitude larger than AM. 
We also investigate the high density behaviour of the equation of state in TM1 model 



Besides scalar self interaction terms, TM1 model includes non-linear u meson term 



The parameters of TM1 model are listed in Table I and Table II. In this model, antikaon 
condensations are found to be second order phase transitions for all values of antikaon optical 
potential depth P Jl9|j20[] . The EoS behaves smoothly even in the high density regime. There 
is no existence of a third family of compact stars in TM1 model because of this smooth EoS. 



V. SUMMARY AND CONCLUSIONS 

We have investigated K~ and A"°-meson condensation in /3-equilibrated hyperonic matter 
within a relativistic mean field model. In this model, baryon-baryon and (anti)kaon-baryon 
interactions are mediated by the exchange of a, uj and p meson. Also, hyperon-hyperon 
interaction is accounted by considering two additional strange mesons a* and (p. In this 
calculation, we have adopted GM1 parameter set. The coupling constants of (anti)kaon- 
meson are determined from a quark model and empirically known values of antikaon optical 
potential depth at normal nuclear matter density. Similarly, the hyperon-meson coupling 
constants are obtained from hypernuclear data and SU(6) symmetry relations. However, the 
coupling constant of hyperons to a* meson is obtained by reproducing the hyperon depth in 
a bath of the hyperons at normal nuclear matter density. 

The condensation of K~ meson in hyperonic matter is described as a first order phase 
transition in this calculation. The /^-equilibrated matter is composed of three phases - pure 
hadronic phase, pure K~ condensed phase and the mixed phase of two pure phases. Along 
with the Gibbs phase rules, global conservation laws are employed in the mixed phase. 

The appearance of antikaon condensation is quite sensitive to the EoS and also depends 
on the value of antikaon optical potential depth. In the presence of hyperons, an EoS becomes 
softer thus delaying the onset of antikaon condensation. For antikaon optical potential depth 
l^x( n o)| < 160 MeV, the early appearance of negatively charged hyperons in particular H~ 
hyperons, reduces the electron chemical potential. And the antikaon condensation does 
not occur even at 7.5n . In our calculation with U^(no) = —160 MeV, the first order 
phase transition to K~ condensation begins at 4.0no- The pure hadronic phase does not 
contain any hyperon whereas only A hyperons appear in the mixed phase. As soon as K~ 
condensation sets in, it is energetically favourable for K~ mesons in the zero momentum 
state to maintain charge neutrality in the system. Consequently, the formation of negatively 
charged hyperons are delayed to higher density. Immediately after the end of the mixed 
phase, K° condensation occurs in the pure K~ condensed phase. Here, we have treated K° 
condensation as a second order phase transition. With the formation of K° condensation, 
abundances of neutron and protons become equal. Negatively charged hyperons - S~ and 
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E~ start populating the system around ~ 7n . And the density of K~ condensate dwindles 
because hyperons having conserved baryon number are energetically more favourable for 
charge neutrality. Consequently, K° condensate dominates the high density regime. 

The equation of state becomes softer in the presence of antikaon condensate compared 
with the situation without any condensate. Among all EoS considered here, the softest 
EoS is the one that includes both K~ and K° condensate. This softening leads to the 
reduction in maximum masses of compact stars. The behaviour of the EoS with both 
antikaon condensates at very high density ~ 8 — 10n and its connection to the possible 
existence of stable super dense stars beyond the neutron star branch have been investigated 
in this work. With the EoS corresponding to hyperonic matter with K~ and K° condensate, 
we find, after the termination of the neutron star branch, an unstable region followed by 
a stable sequence of compact stars called the third family of compact stars. It is noted 
that neutron star twins could exist because of partially overlapping mass regions of the 
neutron star branch and the third family branch. The high density twin in the third family 
branch has a different composition and smaller radius than its low density counter part. The 
observation of two compact stars with almost similar masses but different radii may shed 
light on the existence of a third family of super dense stars. 



13 



REFERENCES 



N.K. Glendenning, Compact Stars (Springer- Verlag, New York, 1997). 
U.H. Gerlach, Physical Review 172, 1325 (1968). 

B.K. Harrison, K.S. Thorne, M. Wakano and J. A. Wheeler, Gravitation Theory and 
Gravitational Collapse, (University of Chicago Press, Chicago, 1965). 
F. Weber, Pulsars as Astrophysical Laboratories for Nuclear and Particle Physics (IOP 
Publishing, Bristol, 1998). 

N.K. Glendenning and C. Kettner, Astron. Astrophys. 353, L9 (2000). 

K. Schertler, C. Greiner, J. Schaffner-Bielich and M.H. Thoma, Nucl. Phys. A677, 463 

(2000). 

J. Schaffner-Bielich, M. Hanauske, H. Stocker and W. Greiner, |astro-ph / 000549C . 



E.S. Fraga, R.D. Pisarski and J. Schaffner-Bielich, D 63, 121702 (2001). 
N.K. Glendenning and J. Schaffner-Bielich, Phys. Rev. C 60, 025803 (1999). 

D. B. Kaplan and A.E. Nelson, Phys. Lett. B 175, 57 (1986); 

A. E. Nelson and D.B. Kaplan, Phys. Lett. B 192, 193 (1987). 

G.E. Brown, K. Kubodera, M. Rho and V. Thorsson, Phys. Lett. B 291, 355 (1992). 

V. Thorsson, M. Prakash and J.M. Lattimer, Nucl. Phys. A572, 693 (1994). 

P.J. Ellis, R. Knorren and M. Prakash, Phys. Lett. B 349, 11 (1995). 

M. Prakash, I. Bombaci, M. Prakash, P.J. Ellis, J.M. Lattimer and R. Knorren, Phys. 

Rep. 280, 1 (1997). 

T. Muto, Prog. Theo. Phys. 89, 415 (1993). 

R. Knorren, M. Prakash and P.J. Ellis, Phys. Rev. C 52, 3470 (1995). 
J. Schaffner and I.N. Mishustin, Phys. Rev. C 53, 1416 (1996). 

J.A. Pons, S. Reddy, P.J. Ellis, M.Prakash, and J.M. Lattimer, Phys. Rev. C 62 035803 
(2000). 

S. Pal, D. Bandyopadhyay and W. Greiner, Nucl. Phys. A674, 553 (2000). 
S. Banik and D. Bandyopadhyay, Phys. Rev. C 63, 035802 (2001). 
J. Boguta and A.R. Bodmer, Nucl. Phys. A292, 413 (1977). 

B. D. Serot and J.D. Walecka, Adv. in Nucl. Phys. 16, 1 (1986). 
N.K. Glendenning, Phys. Rev. D 46, 1274 (1992). 

N.K. Glendenning and S.A. Moszkowski, Phys. Rev. Lett. 67, 2414 (1991). 
R.E. Chrien and C.B. Dover, Annu. Rev. Nucl. Part. Sci. 39, 113 (1989). 

C. B. Dover and A. Gal, Prog. Part. Nucl. Phys 12, 171 (1984). 
T. Fakuda et al, Phys. Rev. C 58, 1306 (1998). 

P. Khaustov et al, Phys. Rev. C 61, 054603 (2000). 
R.S. Hayano et al., Phys. Lett. B 231, 355 (1989). 

E. Friedman, A. Gal and C.J. Batty, Nucl. Phys. A579, 518 (1994); 
C.J. Batty, E. Friedman and A. Gal, Phys. Rep. 287, 385 (1997). 

J. Schaffner, C.B. Dover, A. Gal, D.J. Millener, C. Greiner and H. Stocker, Ann. Phys. 
(N.Y.) 235, 35 (1994). 

E. Friedman, A. Gal, J. Mares and A. Cieply, Phys. Rev. C 60, 024314 (1999). 

V. Koch, Phys. Lett. B 337, 7 (1994). 

T. Waas and W. Weise, Nucl. Phys. A625, 287 (1997). 

G.Q. Li, C.-H. Lee and G.E. Brown, Phys. Rev. Lett. 79, 5214 (1997); Nucl. Phys. 
A625, 372 (1997). 



14 



[36] S. Pal, CM. Ko, Z. Lin and B. Zhang, Phy. Rev. C 62, 061903(R) (2000). 
[37] G. Baym, C.J. Pethick and P. Sutherland, Astrophys. J. 170, 299 (1971). 
[38] J.W. Negele and D. Vautherin, Nucl. Phys. A207, 298 (1974). 

[39] C. Kettner, F. Weber, M.K. Weigel and N.K. Glendenning, Phys. Rev. D 51, 1440 
(1995). 

[40] J.M. Bardeen, K.S. Thorne and D.W. Meltzer, Astrophys. J. 145, 505 (1966). 
[41] D.W. Meltzer and K.S. Thorne, Astrophys. J. 145, 514 (1966). 
[42] Y. Sugahara and H. Toki, Nucl. Phys. A579, 557 (1994). 
[43] A.R. Bodmer, Nucl. Phys. A526, 703 (1991). 



15 



TABLES 



TABLE I. The nucleon-meson coupling constants in the GM1 set are taken from Ref. [24]. In 
this relativistic model, the baryons interact via nonlinear cr-meson and linear w-meson exchanges. 
The coupling constants are obtained by reproducing the nuclear matter properties of binding energy 
E/B = —16.3 MeV, baryon density no = 0.153 fm~ 3 , asymmetry energy coefficient a asy = 32.5 
MeV, incompressibility K = 300 MeV, and effective nucleon mass m* N /m^ = 0.70. The hadronic 
masses are m N = 938 MeV, m a = 550 MeV, = 783 MeV, and m p = 770 MeV. The parameter 
set TM1 is obtained from Ref. [42] which incorporates nonlinear exchanges in both a and co mesons. 
The nuclear matter properties in the TM1 set are E/B = —16.3 MeV, tiq = 0.145 fm -3 , a asy = 36.9 
MeV, K = 281 MeV, and m* N /mN = 0.634. All the hadronic masses in this model are same as 
GM1 except for a-meson which is m a = 511.198 MeV. All the parameters are dimensionless, except 
Q2 which is in fm -1 . 





9aN 


9wN 


9pN 


92 


53 


54 


GM1 


9.5708 


10.5964 


8.1957 


12.2817 


-8.9780 




TM1 


10.0289 


12.6139 


4.6322 


-7.2325 


0.6183 


71.3075 



TABLE II. The coupling constants for antikaons (K) to cr-meson, g a x, for various values of K 
optical potential depths (in MeV) at the saturation density. The results are for the GM1 and 



TM1 set. 
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-140 


-160 


-180 


GM1 


0.9542 


1.6337 


2.3142 


2.9937 


3.6742 


TM1 


0.2537 


0.8384 


1.4241 


2.0098 


2.5955 
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FIG. 1. The proper number densities rii of various compositions in /3-equilibrated hyperonic matter 
including K~ condensate for GM1 model and antikaon optical potential depth at normal nuclear 
matter density = —160 MeV. 
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FIG. 2. The proper number densities of various compositions in /3-equilibrated hyperonic matter 
including both K~ and K° condensate for GM1 model and antikaon optical potential depth at 
normal nuclear matter density = —160 MeV. 
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FIG. 3. The equation of state, pressure P vs. energy density e in GM1 model. The results are 
for hyperonic matter (dashed line), hyperonic matter including K~ condensate (dotted line) and 
K~ and K° condensation in hyperonic matter (solid line) and in hyperonic matter excluding S 
(dashed-dot line) calculated with the antikaon optical potential depth at normal nuclear matter 
density of U R = -160 MeV. 
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FIG. 4. The compact star mass sequences are plotted with central energy density for GM1 model 
and the antikaon optical potential depth of = —160 MeV. The star masses of hyperonic matter 
with K~ condensate and with further inclusion of K° condensate are shown here. In the latter 
C MS C , M second sequence of compact stars appears in the high energy density regime. 
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FIG. 5. The mass-radius relationship for compact star sequences for hyperonic matter with K~ 
condensate and with further inclusion of K° condensate for GM1 model and antikaon optical 
potential depth of Ug = —160 MeV. The mass-radius relationship for the third family branch is 
shown in the inset. 
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